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Abstract 

Taking into account the SU(3) f breaking effects, the strong coupling constants 
of the 7r, K and 77 mesons with decuplet baryons are calculated within light cone 
QCD sum rules method. It is shown that all coupling constants, even in the case of 
SU(3)f breaking, are described in terms of only one universal function. It is shown 
that for H*° — > H*°r/ transition violation of SU(3)f symmetry is very large and for 
other channels when SU(3) t symmetry is violated, its maximum value constitutes 
10% - 15%. 
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1 Introduction 



Excited experimental results are obtained on pion and kaon photo- and electric-production 
of nucleon during last several years. These experiments are performed at different centers, 
such as MAMI, MIT, Bates, BNL and Jefferson Laboratories. To study the properties of 
the resonances from the existing data, the coupling constants of ir, K and rj mesons with 
baryon resonances are needed. 

In extracting the properties of baryon resonances, the hadronic reactions also play an 
important role. Therefore, for a more accurate description of the experimental data, re- 
liable determination of the strong coupling constants of pseudoscalar mesons is needed. 
Calculation of the strong coupling constants of pseudoscalar mesons with baryons (BBP) 
using the fundamental theory of strong interactions, QCD, constitutes a very important 
problem. The strong coupling constants of BBP belongs to the nonperturbative sector of 
QCD and for estimating these couplings, we need som nonperturbative approachs. Among 
all nonperturbative approaches, the most predictive and powerful one is the QCD sum rules 
method [1]. In the present work, we calculate the strong coupling constants of the pseu- 
doscalar mesons with the decuplet baryons within the framework of the light cone QCD 
sum rules (LCSR) method. In this method, the operator product expansion is performed 
over twist rather than dimension of the operators, which is carried out in the traditional 
sum rules. In the LCSR, there appears matrix elements of the non-local operators between 
the vacuum and the corresponding one-particle state, which are defined in terms of the, 
so-called, distribution amplitudes (DAs). These DAs are the main nonperturbative pa- 
rameters of the LCSR method (more about LCSR can be found in [2,3]). Note that, the 
coupling constants of pseudoscalar and vector mesons with octet baryons is investigated 
within the framework of the LCSR in [4], and [5], respectively. 

The paper is organized as follows. In section 2, the the strong coupling constants of the 
pseudoscalar mesons with the decuplet baryons are calculated within the framework of the 
LCSR method, and relations between these coupling constants are obtained where SU (3) j 
symmetry breaking takes place. In section 3, the numerical analysis of the obtained sum 
rules for the pseudoscalar-meson decuplet-baryon coupling constants is performed. 

2 Light cone QCD sum rules for the pseudo scalar- 
meson decuplet-baryon coupling constants 

In this section, we obtain LCSR for the pseudoscalar-meson decuplet-baryon coupling 
constants. For this aim, we consider the following correlation function 

Kr B2V = i J d*xe^(V(q) \T{T,*(x)tf(0)}\0) , (1) 

where V(q) is the pseudoscalar-meson with momentum q, r/^ is the interpolating current 
of the considered decuplet baryon. The sum rules for the above-mentioned correlation 
function can be obtained, on the one side, by calculating it in terms of the physical states 
of hadrons (phenomenological part), and on the other side, calculating it at p 2 — > — oo in 
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the deep Euclidean region in terms of quarks and gluons (theoretical part), and equating 
both representations through the dispersion relations. 

Firstly, let us concentrate on the calculation of the phenomenological side of the cor- 
relation function (1). The phenomenological part can be obtained by inserting a complete 
set of baryon states having the same quantum numbers as the interpolating current rj^. 
Isolating the ground state of baryons, we obtain 

IE— = I + , (2) 

P 2 " L 2 Pi '"l 

where p 1 = p 2 + q, rrii is the mass of baryon Bi, and • • ■ represents the contributions of the 
higher states and the continuum. 

The matrix elements of the interpolating current between vacuum and the hadron states 
is determined as: 

(0 1^1 B(p, s)> = A B u„(p, s) , (3) 

where Xb is the overlap amplitude, and u^p, s) is the Rarita-Schwinger tensor spinor with 
spin s. The matrix element (B 2 (p 2 )V(q) | Bi(pi)) is parametrized as: 



(B 2 (p 2 )V(q) | B 1 ( Pl )) = g Bl B 2 vu a (p2h 5 u a ( Pl ) , (4) 

In order to obtain the expression for the phenomenological part of the correlation func- 
tion, the summation over the spins of the Rarita-Schwinger fields is performed, i.e., 

u„(p, s)u v (p, s) = W + m)(- + ^7„ 7 „ ~ -J^r ~ '^"^ ) ■ ( 5 ) 

s ^ / 

In principle, Eqs.(2-5) allow us to write down phenomenological part of the correlation 
function. However, here following two principal problems appear: 1) not all Lorentz struc- 
tures are independent; 2) not only spin-3/2, but also spin-1/2 states contribute. Indeed, 
the matrix element of the current 77^, sandwiched between the vacuum and the spin-1/2 
states, is different than zero and determined in the following way: 

(0 |^| B(p, s = 1/2)) = A{Ap, - m llt )u{p, s = 1/2) , (6) 

where the condition ^^rf = has been used. 

There are two different alternatives to remove the unwanted spin-/ 1/2 contribution 
and take into account only the independent structures: 1) ordering the Dirac matrices in 
a specific way and eliminate the ones that receive contributions from spin-/ 1/2 states; 2) 
introduce projection operators for the spin-3/2, that do not contain spin-1/2 contribution. 

In the present work, we have used the first approach and choose the ^^1^15 ordering 
of the Dirac matrices. Having chosen this ordering for the Dirac matrices, we obtain 

Iin U = 5 — ^ B ^ BiV „ (guvp'ilh + other structures with 7„ at the beginning and 
[pi - m{){pi - mi) V 

7j, at the end, or terms that are proportional to p\ v or p 2 A . (7) 
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The advantage of choosing the structure g^vfil^ is in the fact that, the spin-1/2 states 
do not give contribution to this structure. This fact immediately follows from Eq. (6), 
which tells that spin-1/2 states contribution is proportional to or 7^. 

In order to calculate the theoretical part of the correlation function (1) from the QCD 
side, we need the explicit expressions of the interpolating currents of the decuplet baryons. 
The interpolating currents have the following forms [6] : 



abc 



?r indices and C is the charge conjugation open 



(8) 



where a, b, c are the color indices and C is the charge conjugation operator. The values of 
A and the quark flavors q±, q 2 and g 3 for each decuplet baryon are presented in Table-1. 
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Table 1: The values of A and the quark flavors q±, q 2 and q 3 

Before presenting detailed calculation of the correlation function from the QCD side for 
determination of the coupling constants of pseudoscalar-mesons with decuplet-baryons, let 
us establish the relation among the correlation functions, more precisely, relations among 
the coefficients of the invariant functions for the structure g^ilh- F° r this aim, we will 
follow the works of [4, 7] , and we will show that all correlation functions which describe the 
strong coupling constants of pseudoscalar-mesons with decuplet-baryons can be written in 
terms of only one invariant function. It should especially be noted that the approach we 
present below automatically takes into account the SU (3) / symmetry breaking effects. 

In obtaining relations among the invariant functions, similar to works [4,5], we start by 
considering the the correlation function describing £*° — > S*°7r° transition, this correlation 
function can formally be written in the following form: 

n E *°^ s * 07r ° = j„ni(n, d, s) + w<w ni(u, d, s) + g nss U 2 (u : d, s) , (9) 
where 7r° current can formally be written as 

J = 9n qq qihq , (10) 

q=u,d,s 
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where g w o uu = —g^dd — I/a/2 and <7„-o ss = for ir° meson. The functions III, 11^ and II2 
describe radiation 7r° meson from u, d and s quarks of the S*° baryon, respectively. 

The interpolating current 77 s * is symmetric under the change u <-> d, and therefore 
H'i(u, d, s) = H\(d, u, s). Hence, Eq. (9) can be written as: 



*0 ,v*0_0 1 



-j-|-£*u_+X;* u 7r' 



Ili(u,d,s) - ni(d,u,s)] . (11) 



V2 

For convenience, let us introduce the notations 

ni(u,d,s) = (m« 0) , 

n 2 (u,d,s) = <ss|S*°S*°|0) . (12) 

Obviously, Il 2 = for the transition S*° -> £*°7r°. 

In the transition with 77 meson, the situation is more complicated, since strange quark 
is in the quark content of the rj meson. In the present work, we neglect the mixing between 
the rj and rj' mesons and 77 meson current is taken to have the following form: 

J v = -^=(m7 5 m + d"/ 5 d - 2s7 5 s) . (13) 

A simple analysis shows that the S*° — > transition has the similar form as is given in 
Eq. (9) 

n E *°^ E * 0?J = g^Jl^u, d, s) + g r)dd U' l {u, d, s) + g vss U 2 (u, d, s) . (14) 
Using the definition given in Eq. (12), one can easily show that 

H 2 (u,d, s) = Ui(s,d,u) . (15) 
For this reason, using Eqs. (13) and (15), we get from Eq. (14) 

n E^ E .o„ = 1 r n ^ Uj dj s ^ + n ^ Uj s ^ _ 2n ^ Sj dj J _ ^ 

V6 L -I 

The invariant function describing the S* + — > S* + 7r° transition can be obtained from Eq. 
(9) with the help of the replacements d — > u in IIi(u, d, s) and using the fact £*° = -^2*+, 
which results in 

4n(u,u,s) = 2{uu | £*+£*+ 1 0) . (17) 

The presence of factor 4 on the left-hand side of Eq. (17) can be explained as follows. 
Each S* + contains two u quarks and therefore there are 4 ways that the ir° meson can be 
radiated. Since £*+ does not contain d quark, for the S*° — > S*°7r° transition, it can be 
written from Eq. (9) that 

= v^ni(u,«,s) . (18) 
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The result for the X* — > S* n° transition can easily be obtained by making the re- 
placement u — > d in Eq. (9) and using X*°(w — > d) — v^ffi* - , from which we obtain 

n s '^ s, "'° = g n o dd (dd\Z*-Z*-\0) + g7T o ss (ss\Z*-Z*-\0) 

= -y/mi{d,d,s) . (19) 

In the case of exact isospin symmetry, it follows from Eqs. (11), (18) and (19) that 

nS *o_ >s *o 7r o = u x« + ^ +7V o = _ nS *-^ s - 7r o_ 

Let us now calculate the invariant function responsible for the A + — > A + 7r° transition. 
Since A+ = S*+(s -> d), we get from Eq. (18) 

n A+ ^ A+w ° = ^o uu (uu | £*+£*+ 1 0) (s - d) + ^ M (ss | £*+£*+ 1 0) (s - d) 

(u,u,d) ^=rii(d, u) . (20) 

v2 

Similarly, it is not difficult to obtain the relations for the transitions in which A , A ++ 
and A~ decuplet baryons and 7r° meson participate: 



-V / 21I 1 ((i, d, -u) + -^=rii(u, d, d) , 
V2 



n A++ - A++ -° = n s * + -^°( s - u) 

3 

= -^IIi(«, u,u) , 

n A -- A -« = n 53 -^-*°( a -d) 

3 

= --^=IIi(d,cZ,d) , 

n H ^ H ^° = --J=n 1 (d, s , s ). (21) 

We can proceed now to obtain similar relations in the presence of charged n meson. 
In order to obtain these relations, we consider the matrix element (dd 0), where d 

quarks from each S*° form the final dd state and, u and s quarks are the spectators. In 
the matrix element {ud |E* + E*°| 0), d quark from S*° and u quark from S* + form the ud 
state and the other u and s quarks are the spectators. For these reasons, it is natural to 
expect that these matrix elements should be proportional to each other. Direct calculations 
confirm this expectation, i.e., 

n2 *o^ +7r - = («d|£*+£*°|o) = v^(dd|£*°£*°|0) 

= V2U[ (u, d, s) = V^U^d, u,s) . (22) 
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Making the replacement u <-> d in Eq. (22), we get 

= y/m 1 {u,d,s) . (23) 

Along the same lines of reasoning, similar calculations for for A and S decuplet baryons 
are summarized below: 

-V2(fm|S*°S*°|0) = ni(d, s,s) , 

rii(u, s,s) , 





= (du\E*°E*-\0) 
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J-JA+-+A + + 7T- 


= v / 3II 1 (-u, -u, u) , 


j-jA-^A°7r- 


= >/3ni(d,d,d) . 



(24) 

The correlation function involving K meson can be obtained from the previous results 
as follows: 



U s *-^*- K ° = U-°^ +K -(u - d) = 2U 1 (d,d,s) 

u x* + ^~*o K+ = uS *o^ +K - (M ^ s) = 2Hl ( S) s? M ) . (25) 

Remaining correlation functions involving ir and K mesons are presented in Appendix-A. It 
follows from the results presented above that, all coupling constants of pseudoscalar mesons 
with decuplet baryons can be expressed by only one independent invariant function, which 
constitutes the main result of the present work. 

Having obtained this result, our next task is the calculation of the correlation function 
from QCD side. The correlation function in deep Euclidean domain p\ — > — oo, p\ — > — oo, 
can be calculated using the operator product expansion. For this purpose the propagators 
of light quarks, as well as their distribution amplitudes (DA's) are needed. The matrix 
elements (V(q) \q(xi)Tq'(x2)\ 0) that parametrized in terms of DA's are given in [8-10]: 

(V(p)\q(x) W q(0)\0) 

(V(p)\q(x)i l5 q(0)\0) 
(V(p)\q(x)a af3l5 q(0)\0) 
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Jo 
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(26) 
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and qi and g 2 are the quarks in the meson V, Va = dagda g da g S(l — a q — a q — a g ), and 
and the DA's ip-p(u), A(u), B(-u), ip P (u), (p a (u), T(aij), A±(cti), A\\(cxi), V±(aii) and Vy(ai) 
are functions of definite twist and their expressions are given in the next section. 

For the calculation of the correlation function, we use the following expression for the 
light quark propagator, 
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(27) 



where 7# ~ 0.577 is the Euler Constant. In the numerical calculations the scale parameter 
A is chosen as factorization scale, i.e., A = 0.5 4- 1.0 GeV. This point is discussed in detail 
in [11,12]. 

Using Eqs. (26) and (27) and separating the coefficient of the structure g^ilb, the 
theoretical part of the correlation function can be calculated straightforwardly. Equating 
the coefficients of the structure g^ilb from physical and theoretical parts, and performing 
Borel transformation in the variables p\ = p 2 and p\ = (p + q) 2 in order to suppress the 
higher states and continuum contributions [13, 14], we get the sum rules for the correspond- 
ing pseudoscalar-meson decuplet-baryon coupling constants. 

As the result of our calculations, we obtain the following expression for the invariant 
function IIi(-u, d, s). 
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(28) 



In calculating the coupling constants of pseudoscalar-mesons with decuplet baryons, 
the value of the overlap amplitude of the hadron is needed. This overlap amplitude is 
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determined from the analysis of the two-point function which is calculated in [13, 14]. Our 
earlier considerations reveal that the interpolating currents of decuplet baryons can all be 
obtained from the S*° current, and for this reason we shall present the result only for the 
overlap amplitude of £*°: 
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where x = s /M 2 . 

The contribution of the higher states and continuum in ITx are subtracted by taking 
into account the following replacements 
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(30) 



3 Numerical analysis 

In this section, we present the numerical calculations for the sum rules for the couplings of 
the pseudoscalar-mesons with decuplet-baryons. The main nonperturbative parameters of 
LCSR are the DA's of the pseudoscalar mesons, whose explicit forms entering Eq. (26) are 
given in [8-10]: 
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(32) 



The values of the parameters a[, a 2 , 773, 774, u>3, and W4 entering Eqs. (32) are given in 
Table (2) for the ir, K and rj mesons. 

In the numerical calculations, we set Mf = Mf = 2M 2 due to the fact that the masses of 
the initial and final baryons are close to each other. With this choice, we have Mo = 1/2. The 
values of the other input parameters entering the sum rules are: (qq) = — (0.24±0.01 GeV) 3 , 



(0.8 ± 0.2) GeV 2 [13], f n = 0.131 GeV, f K = 0.16 GeV, and /„ = 0.13 GeV 
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Table 2: Parameters of the wave function calculated at the renormalization scale [i — 1 GeV 

The sun rules for the coupling constant of pseudoscalar-mesons with decuplet-baryons 
contain two auxiliary, namely, Borel parameter M 2 and the continuum threshold sq. Ob- 
viously, we need to find such regions of these parameters where coupling constants are 
practically independent of them. 

The upper limit of M 2 can be found by requiring that the higher states and continuum 
contributions to the correlation function should be less than 40-50% of the total value of 
the correlation function. The lower bound of M 2 can be obtained by demanding that the 
contribution of the highest term with power 1/M 2 is less than, say, 20-25% of the highest 
power of M 2 . Using these two conditions, one can find regions of M 2 where the results for 
the coupling constants are insensitive to the variation of M 2 . The value of the continuum 
threshold is varied between 2.5 GeV 2 < s < A GeV 2 . 

As an example, in Fig. (1), we depict the dependence of the g s * + ^ s * +7r ° coupling 
constant on M 2 at fixed values of the continuum threshold. From this figure, one can see 
that the g E * + ^ s * +7r ° coupling constant demonstrates good stability to the variation in M 2 . 
The numerical results for the coupling constants of pseudoscalar-mesons with decuplet- 
baryons are presented in Table (3). Note that in this Table, we give only those results 
which are not obtained from each other by SU (2) and isotopic spin relations. It should be 
remembered that the sum rules cannot fix the signs of the residues and for this reason the 
signs of the couplings are not fixed. However, they can be fixed if we use SU (3) / symmetry 
(more about this issue, see [4]). From this Table, we can deduce the following conclusions: 

• In all considered couplings except S*° — > T,* !] our predictions consist with the SU (3) / 
symmetry. Maximum violation of SU (3) / symmetry is about 15%. 

• In SU(3)f symmetry limit coupling constant for S*° — > S*°r] transition is equal to 
zero, but our prediction on this constant differs from zero considrably when violation 
of £77(3) f symmetry is taken into account. Only for this channel violation of SU (3) j 
symmetry is huge. In principle, investigation of this coupling constant can shed light 
on the structure of rj meson. 

• Sign of coupling constant of decuplet baryons to K meson and also S*° — > E*°r] is 
negative, but for all other cases is positive. 

In summary, cosidering the SU(3) f symmetry breraking effects, the coupling constants of 
the decuplet baryons with pseudoscalar n, K and rj mesons have been calculated in the 
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framework of light cone QCD sum rules. It was shown that all aforementiond coupling 
constants is described with the help of one universal function. We obtained that for S*° — > 
S*°77 transition, violation of SU (3) / is very large. 



channel 


Coupling 


Coupling in SU(3) limit 




11.3±1.2 


11.0±1.2 


A+^A+^° 


5.5±0.8 


5.5±0.8 


>S* 7T 


5.2±0.3 


5.5±0.5 




-17.4±1.6 


-18.0±1.8 


■z:*0 y Yj*~^~ K~ 


-25.4±2.5 


-27.0±2.8 


S*+^A++K- 


-21.2±1.8 


-22.0±2.0 




-20.7±1.6 


-22.2±2.0 




-22.2±1.6 


-27.0±2.8 




0.65±0.05 


O.OiO.O 




12.5±1.5 


12.6±1.5 




-11.2±1.0 


-13.2±1.8 



Table 3: Coupling constants of pseudoscalar-mesons with decuplet-baryons 
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Appendix A 



In this appendix we present the correlation functions involving ir, K and r] mesons which 
is not given in the main text. 

• Correlation functions for the couplings involving n + meson 

• Correlation functions for the couplings involving tt~ meson 

n A °^ A+7r " = 2U 1 (u,u,d) . 

• Correlation functions for the couplings involving K meson 

n £*^A+K- =v /2 ni ( S ,M), 

n E*+^A++^- = ^(u^u) , 
n E*-^A°tf- =ni ( S ,d,d) , 

n n-^"K- = ^3 ni ( S ,s,s) , 

n AO^E* + ^- = Q ^ 

n A°^E*"* + =ni ( a ,d,d) , 
nS * + _ >A o x+ = Q ^ 

^0^,0^0 = ^.o^ojfo = nS .o^o = nS .o_ H ,o / ,o = ^ ^ ^ 

n E-->H*-i?o = n --^E-^ = 2Ui{s s d) ) 

n n-^s-ff° = uE *-^ n - R o = uE ^ n - K+ = u n-^-K" = jjB-^n-^ = V3 ni ( s , s , s ) , 

nS * + ^ A+ X0 = n A + ^ + X0 = ni(aj1i>u) j 
n A-^-X« = nS *-^ A -^o = ^ 

n A + ^E* + *° = n E* + ^A + ^ = ni ( S , M , M ) , 
n A-^E*-i.o = U ^ A - K0 = ^ Ui{s d d) j 
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• Correlation functions for the couplings involving 7] meson 



n A°-+A°r; 



s .o_ >3 .o^ 



1 

71 

2 

^6 
2 

1 

^6 
>/6 



Il^-u, d, s) + Ui(d, u, s) — 2n!(s, d, u)J , 

IIi(u,U,s) - Il 1 (s,M,M)J , 

d, s) — d, d)J , 
2Ui(u, u, d) + IIi (<i, u, u) 



rr 



\/6 

^-n!(d,d,d) , 
i 

71 
i 

1 

V6 



2IIi(d, d, u) + Hi(u, d, d) 
IIi(u,s,s) -4IIi(s, s,u)] , 
n^d, s, s) - 4IIi(s, s, d)j . 
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Figure 1: The dependence of the g s * + ^ s * +7r ° coupling constant on M 2 at fixed values of 
the continuum threshold. 
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